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Abstract 

In the era of high precision CMB measurements, systematic effects are beginning to 
limit the ability to extract subtler cosmological information. The non-circularity of 
the experimental beam has become progressively important as CMB experiments 
strive to attain higher angular resolution and sensitivity. The effect of non-circular 
beam on the power spectrum is important at multipoles larger than the beam-width. 
For recent experiments with high angular resolution, optimal methods of power spec- 
trum estimation are computationally prohibitive and sub-optimal approaches, such 
as the Pseudo-C/ method, are used. We provide an analytic framework for correct- 
ing the power spectrum for the effect of beam non-circularity and non-uniform sky 
coverage (including incomplete/masked sky maps). The approach is perturbative 
in the distortion of the beam from non-circularity allowing for rapid computations 
when the beam is mildly non-circular. When non-circular beam effect is important, 
we advocate that it is computationally advantageous to employ 'soft' azimuthally 
apodized masks whose spherical harmonic transform die down fast with m. 
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1 Introduction 

The fluctuations in the Cosmic Microwave background (CMB) radiation are 
theoretically very well understood, allowing precise and unambiguous predic- 
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tions for a given cosmological model (Jig |2j). Consequently, measurement of 
CMB anisotropy has spearheaded the remarkable transition of cosmology into 
a precision science. The transition has also seen the emergence of data analy- 
sis of large complex data sets as an important and challenging component of 
research in cosmology. Increasingly sensitive, high resolution, measurements 
over large regions of the sky pose a stiff challenge for current analysis tech- 
niques to realize the full potential of precise determination of cosmological 
parameters. The analysis techniques must not only be computationally fast to 
contend with the huge size of the data, but, the higher sensitivity also limits 
the simplifying assumptions that can be then invoked to achieve the desired 
speed without compromising the final precision goals. There is a worldwide 
effort to push the boundary of this inherent compromise faced by the current 
CMB experiments that measure the anisotropy in the CMB temperature and 
its polarization. 

Accurate estimation of the angular power spectrum, Ci, is arguably the fore- 
most concern of most CMB experiments. The extensive literature on this topic 
has been summarized in literature (j^S)- F° r Gaussian, statistically isotropic 
CMB sky, the C\ that corresponds to the covariance that maximizes the mul- 
tivariate Gaussian PDF of the temperature map, AT(q) is the Maximum 
Likelihood (ML) solution. Different ML estimators have beenproposed and 
implemented on CMB data of small and modest sizes 

(US 13 3 Hi), while 

it is desirable to use optimal estimators of C\ that obtain (or iterate toward) 
the ML solution for the given data, these methods are usually limited by 
the computational expense of matrix inversion that scales as iVj with data 
size Nd ([lOt llll ) . Various strategies for speeding up ML estimation have been 
proposed, such as, exploitin g th e symmetries of the scan strategy fl2Th using 
hierarchical decomposition (|13T). iterative multi-grid method etc. Vari- 
ants employing linear combinations of AT(q) such as ai m on set of rin gs i n 
the sky can alleviate the computational demands in special cases |l5[ lid). 
Other promising 'exact' power estimation methods have been recently pro- 
posed JII El G3) 



However there also exist computationally rapid, sub-optimal estimators of C\. 
Exploiting the fast spherical harmonic transform (~ N^ 2 ), it is possible to 
estimate the angular power spectrum C\ = J2 m I a/rJ 2 / (2/ + 1) rapidly (2(1 21). 
This is commonly referred to as the Pseudo-Q method (22). 1 It has been 
recently argued that the need for optimal estimators may have been over- 
emphasized since they are computationally prohibitive at large / . Sub-optimal 
estimators are computationally tractable and tend to be nearly optimal in 
the relevant high I regime. Moreover, already the data size of the current 
sensitive, high resolution, 'full sky' CMB experiments such as WMAP have 



1 Analogous approach employing: fast estimation of the correlation function C(q-q') 
have also been explored (23; 24). 
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been compelled to use sub-optimal Pseudo-C; related methods On the 

other hand, optimal ML estimators can readily incorporate and account for 
various systematic effects, such as noise correlations, non-uniform sky coverage 
and beam asymmetries. The systematic correction to the Pseudo-C; power 
spectrum estimate arising from non-uniform sky coverage has been studied and 
implemented for CMB temperature (27) and polarization (28). The systematic 
correction for non circular beam has been studied by us (|29f). Here we extend 
the results to include non-uniform sky coverage. 



It has been usual in CMB data analysis to assume the experimental beam re- 
sponse to be circularly symmetric around the pointing direction. However, any 
real beam response function has deviations from circular symmetry. Even the 
main lobes of the beam response of experiments are generically non-circular 
(non-axisymmetric) since detectors have to be placed off-axis on the focal 
plane. (Side lobes and stray light contamination add to the breakdown of this 
assumption). For highly sensitive experiments, the systematic errors arising 
from the beam non-circularity become progressively more important. Drop- 
ping the circular beam assumption leads to major complications at every stage 
of analysis pipeline. The extent to which the non-circularity affects the step 
of going from the time-stream data to sky map is very sensitive to the scan- 
strategy. The beam now has an orientation with respect to the scan path that 
can potentially vary along the path. This implies that the beam function is 
inherently time dependent and difficult to deconvolve. 



Even after a sky map is made, the non- circularity of the effective beam affects 
the estimation of the angular power spectrum, Ci, by coupling the power at 
different multipoles, typically, on scales beyond the inverse angular beam- 
width. Mild deviations from circularity can be addressed by a perturbation 



approach (|30t l3l[ ) and the effect of non-circularity on the estimation of CMB 



power spectrum can be studied (semi) analytically (29). Fig. 1 shows the 
predicted level of non-circular beam correction in our formalism for elliptical 
beams with fwhm beam-width of 0.22° compared to the non-circular beam 
corrections computed in the recent data release by WMAP (j4ll ). 



To avoid contamination of the primordial CMB signal by Galactic emission, 
the region adjoining the Galactic plane is masked from maps. If the Galactic 
cut is small enough, then the coupling matrix will be invertible, and the two- 
point correlation function can be determined on all angular scales from the 
data within the uncut sky (32). Hivon et al. ([271 ) present a technique (MAS- 
TER) for fast computation of the power spectrum taking accounting for the 
galactic cut, but for circular beams. In our present work, we present analytical 
expressions for the bias matrix of the Pseudo-Ci estimator for the incomplete 
sky coverage, using a non- circular beam. 
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Fig. 1. The predicted non-circular beam correction for a CMB experiment with ellip- 
tical Gaussian beam with fwhm beam-width of 0.22° and two values of eccentricity, 
e = 0.6 and e = 0.4 are shown as dashed-dot lines (labeled Th. Pr.). The solid 
curves are the non circular beam corrections estimated by the WMAP team for the 
Q,V and W channels. The fwhm beam-width of 0.22° corresponds closely to the 
cleanest and highest resolution W-band beam-width. For larger beam-widths, the 
theoretical curves would roughly shift left in multipole by the ratio of beam-widths. 

2 CMB angular Power spectrum estimation: 

Non-circular beam & Non-uniform coverage effects 

The observed CMB temperature fluctuation is convolved with a beam function 
and contaminated by noise. Further, the CMB signal cannot be obtained for 
full sky because of galactic contamination (and extragalactic point sources). 
We derive the general form of the bias matrix including non-uniform/incomplete 
sky coverage and a general beam function. 

The observed temperature fluctuation field AT(q) is the convolution of the 
"beam" profile _B(q, q') with the real temperature fluctuation field AT(q) 
(ignoring the additive noise term for simplicity) : 



The two point correlation function for a statistically isotropic CMB anisotropy 
signal is 



AT(q) 





C(q,q') = (AT(q)AT(q')> = £ 



(21 + 1) 



Q w,(q> q') 



(2) 



1=0 



where Ci is the angular spectrum of CMB anisotropy signal and the window 



function 



= jd^j (Kit B(qx, q)S(q 2) q')P { (q ■ q'), 



(3) 



encodes the effect of finite resolution through the beam function. A CMB 
anisotropy experiment probes a range of angular scales characterized by a 
window function Wi(q, q'). The window depends both on the scanning strategy 
as well as the angular resolution and response of the experiment. However, 
it is neater to logically separate these two effects by expressing the window 
Wi(q,q!) as a sum of 'elementary' window function of the CMB anisotropy 
at each point of the map (j3o[). For a given scanning strategy, the results can 
be readily generalized using the representation of the window function as sum 
over elementary window functions (see, e.g., (jiot liih). 

For some experiments, the beam function may be assumed to be circularly 
symmetric about the pointing direction, i.e., -B(q, q') = -B(q • q') without 
significantly affecting the results of the analysis. In any case, this assumption 
allows a great simplification since the beam function can then be represented 
by an expansion in Legendre polynomials as 

-| oo 

5 (q-q) = T-E^ + ^^q-q)- (4) 

47r 1=0 

Consequently, it is straightforward to derive the well known simple expression 
W,(q,q') = Bfi^q-q'), (5) 

for a circularly symmetric beam function. 
We define the Pseudo-Q estimator as 

Ci = ±J d% J dOfc tf(qa) U(q 2 ) P,(qi • q 2 ) AT(qi) AT(q 2 ), (6) 

4-7T 47T 



where U(q) denotes the mask function representing the incomplete sky. The 
expectation value of the Pseudo-C/ estimator can be shown to take the form 



(Ci) 



21 



l'=0 



lm=-V 



taqf/(q) Y ln (q) 



Itt 



^y; m (q')B(q,q') 



(7) 



An 
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The integral in the square bracket can be simplified to (|3 



2/ ' 1 E ^Aw< m ,(q,p(q)). (8) 



In 



Att 



m'=-V 



The Beam Distortion Parameter (BDP) (5\ m = 6/ m /6;o is expressed in terms of 



k m = J dO^;(q)S(z,q), 

1 

B, = | ^q-qOP^q-qO^q-q) 



(9) 



where £>(q • q') is the circularized beam obtained by averaging B(z, q) over 
azimuth (/>. Hence, 



5, = / sin 9d9 



Ait 



21 + 1 



^o(q) 



1 2?r 

— / #£?(z,q) 



4tt 



2/ + 1 



fc/o- (10) 



Making a spherical harmonic expansion of the mask function U(q) 

oo I 

U ($) = E E ^lm^m(q) 
i=0 m=-Z 



we can simplify Eq.7 as 

(Q) = E^/'CV- 
v 



The general form of the bias matrix, Aw is thus given by 
Bf {21' + V 1 1 



(11) 



(12) 



4tt (2Z + 



/ n=—l m=—l' 



V oo i" 

E A'm' E E ^"m" 
m'=-Z' !"=0m"=-(" 



ll"l> 



(13) 



where 



jii"i> 

nm"mm' 



= J (Kl^Y ln (q) Y vtmt >(q) £w(q, p(q))- 

47T 



(14) 



To proceed further analytically, we need a model for p(q). We shall continue 
assuming non-rotating beams, i.e. p(q) = 0. We evaluate the integral </" m » mm /, 
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with two different approaches. In the first method, using only the sinusoidal 
expansion of Wigner-<i, we get 



I 11 " 1 ' —?TrS 

nm" mm' 



m" (m—n) " 



I" 



(2Z + l)(2i" + l) ' 
47r kT^-i 



X ^m"M" ( ) ^M"0 



M"=-l" 

n+m+m 1 +m" 



71 



L M0 



/ j a mM' I n I u M'm' 
M'=-V VZ/ 



71 

-1) M+M " +M ' | sin 0d0 e ^ M + M '+ M ") e 



• (15) 



In the alternative method using Clebsch Gordon coefficients, we can evaluate 



~rll I q . 

° nm" mm' ' 



j: 



11" v 

nm"mm' 



;-ir 

Z+Z" 



-m" r 

Cm" 



m" (m—n) 



(2/+ 1)(2Z" + 1) 



x L wo L/ i 
l=|;-/"| 

Z/ 

L' 



L(n+m") 
lnl"m" 



L+V 



E^L' (m-n-m") y-yL'm' 
° L(-n-m")l'm ^LOl'm* 



\L-l'\ 



X 



N= 



V r/ L ' 1-1 d L ' - 

/ , a 0N I r, I "ATm' 1 r> 



\N 



sin 0c?0 e 



IN6 



(16) 



The analytic expressions reduce to the known analytical results for circular 
beam and non-uniform sky covera ge s tudied in ref. (|27l ) and our earlier results 
for non-circular beam for full sky (|29l). 

These results offer the possibility of rapid estimation of the non-circular beam 
effect in the Pseudo-Q estimation. The expression in terms of the d l mm ,(n/2) 
coefficients is the computationally superior approach. These coefficients can 
be computed using stable recurrence relations ( 34t 1351) . In a more detailed 
publication we describe the algorithm in more detail (|39l). The expressions also 
highlight the two aspects to speeding up the computation of the systematic 
effect: 



i. Mildly non-circular beams where the beam distortion parameters (BDP), 
(3i m at each I fall off rapidly with m. This allows us to neglect j3i m for m > 
m beam- For most real beams, mb eam ~ 4 is a sufficiently good approximation. 
This cuts-off the summations over BDP in the expressions for Au>. 

ii. Soft, azimuthally apodized, masks where the coefficients Ui m are small be- 
yond m > m mas k. Moreover, it is useful to smooth the mask in /, such the 
Ui m die off rapidly for I > Z mas k too. 

The mild-circularity perturbation approach has been introduced and discussed 
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Fig. 2. A softened (azimuthally smoothed) mask reconstructed from the WMAP 
Kp2 mask which would reduce the computational cost of estimating the non-circular 
beam effect on the angular power spectrum. 

in ref. (Jsfih . The circularity of the beam has to be addressed in the design 
of the CMB experiments. Our results suggest the systematics due to non- 
circular distortions of the beam are manageable if one ensures the large BDP 
are limited to a few m (i.e., narrow band limited violation of axis-symmetry). 
The beams for many experiments, such as P yth on- V and WMAP are well 
approximated as elliptical Gaussian functions ()3cl 33L 29^ . For radically non- 
axisymmetric beams, modeling the beam in terms of superposition of displaced 
circular Gaussian beams has been proposed (|36t 1371). Our approach allows a 
simple, cost effective extension to modeling with the more general elliptical 
Gaussian beams, or other mildly non-circular beam forms. 

The mask of the galactic region U(q) can be chosen at the time of data 
analysis. The coupling of BDP with Ui m suggests that a judicious choice of 
mask reduces the computational costs of non-circular beam corrections. Fig. 2 
shows a softened version of the Kp2 mask used by the WMAP team (Jjjh , where 
the mask is azimuthally smoothed. The final apodized mask is obtained by 
multiplying an azimuthally smoothed mask raised to a sufficiently large power 
with the original mask and has reduced power at large m (i.e., Ui m is negligible 
for m > m mas k) 2 . In a forthcoming publication we describe the method of 
making soft masks (jUl). For mildly non-circular, nearly azimuthally symmetric 
case, the required number of computation cycle to compute the bias matrix 
up to a multipole I scales as ~ (2mb ea m + l)(2m mas k + 1) / 5 up to leading order 
for I ^> / mas k. Here, m^, CSuia is the cut-off in the beam distortion parameters 



2 Recently, apodized masks (using circularly symmetric smoothing kernels) have 
been recommended in the context of CMB polarization maps (j40|). In our context, 
the mask is predominantly azimuthally smoothed to retain more CMB sky. 



8 



(BDP), f3i m and m mask is the cut-off in Ui m . 



3 Discussion & Conclusion 

The assumptions of non-circular beam leads to major complications at every 
stage of the data analysis pipeline. The extent to which the non-circularity 
affects the step of going from the time-stream data to sky map is very sensitive 
to the scan-strategy. The beam now has an orientation with respect to the scan 
path that can potentially vary along the path. This implies that the beam 
function is inherently time dependent and difficult to deconvolve. 

We extend our analytic approach for addressing the effect of non-circular ex- 
perimental beam function in the estimation of the angular power spectrum C\ 
of CMB anisotropy, which also includes the effect of the galactic cut in the 
entire sky map. Non-circular beam effects can be modeled into the covariance 
functions in approaches related to maximum likelihood estimation 0) and 
can also be included in the Harmonic ring (jlih and ring-torus estimators (["ill). 
However, all these methods are computationally prohibitive for high resolu- 
tion maps and, at present, the computationally economical approach of using 
a Pseudo-C/ estimator appears to be a viable option for extracting the power 
spectrum at high multipoles (J3(). The Pseudo-Q estimates have to be corrected 
for the systematic biases. While considerable attention has been devoted to 
the effects of incomplete/non-uniform sky coverage, no comprehensive or sys- 
tematic approach is available for non-circular beam. The high sensitivity, 'full' 
(large) sky observation from space (long duration balloon) missions have alle- 
viated the effect of incomplete sky coverage and other systematic effects such 
as the one we consider here have gained more significance. Non-uniform cov- 
erage, in particular, the galactic masks affect only CMB power estimation at 
the low multipoles. The analysis accompanying the recent second data from 
WMAP uses the hybrid strategy (0) where the power spectrum at low multi- 
poles is estimated using optimal Maximum Likelihood methods and Pseudo-Q 
are used for large multipoles (fill li^ ). 

The non-circular beam is an effect that dominates at large I comparable to the 
inverse beam width fj29h . For high resolution experiments, the optimal max- 
imum likelihood methods which can account for non-circular beam functions 
are computationally prohibitive. In implementing the Pseudo-C; estimation, 
we have included both the non-circular beam effect and the effect of non- 
uniform sky coverage. Our work provides a convenient approach for estimating 
the magnitude of these effects in terms of the leading order deviations from a 
circular beam and azimuthally symmetric mask. The perturbation approach is 
very efficient. For most CMB experiments the leading few orders capture most 
of the effect of beam non- circularity (Jicft. Our results highlight the advantage 
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of azimuthally smoothed masks (mild deviations from azimuthal symmetry) in 
reducing computational costs. The numerical implementation of our method 
can readily accommodate the case when pixels are revisited by the beam with 
different orientations. Evaluating the realistic bias and error-covariance for a 
specific CMB experiment with non-circular beams would require numerical 
evaluation of the general expressions for A\y using real scan strategy and ac- 
count for inhomogeneous noise and sky coverage, the latter part of which has 
been addressed in this present work. 



It is worthwhile to note in passing that that the angular power C\ contains 
all the information of Gaussian CMB anisotropy only under the assumption 
of statistical isotropy. Gaussian CMB anisotropy map measured with a non- 
circular beam corresponds to an underlying correlation function that violates 
statistical isotropy. In this case, the extra information pres ent ma y b e mea- 
surable using, for example, the bipolar power spectrum Even 
when the beam is circular the scanning pattern itself is expected to cause a 
breakdown of statistical isotropy of the measured CMB anisotropy (|27l). For a 
non-circular beam, this effect could be much more pronounced and, perhaps, 
presents an interesting avenue of future study. 



In addition to temperature fluctuations, the CMB photons coming from dif- 
ferent directions have a random, linear polarization. The polarization of CMB 
can be decomposed into E part with even parity and B part with odd par- 
ity. Besides the angular spectrum Cf T , the CMB polarization provides three 
additional spectra, Cf E , Cf E and Cf B which are invariant under parity trans- 
formations. The level of polarization of the CMB being about a tenth of the 
temperature fluctuation, it is only very recently that the angular power spec- 
trum of CMB polarization field has been detected. The Degree Angular Scale 
Interferometer (DASI) has measured the CMB polarization spectrum over lim- 
ited band of angular scales in late 2002 (|47l). The DASI experiment recently 
published 3-year results of much refined measurements (|48h. More recently, the 
BOOMERanG collaboration reported new measurements of CMB anisotropy 
and polarization spectra (ji^). The WMAP mission has also measured CMB 
polarization spectra (|50t 1511 1 . Correcting for the systematic effects of a non- 
circular beam for the polarization spectra is expected to become important. 
Extending this work to the case CMB polarization is another line of activity 
we plan to undertake in the near future. 



In summary, we have presented a perturbation framework to compute the 
effect of non-circular beam function on the estimation of power spectrum of 
CMB anisotropy taking into account the effect of a non-uniform sky cover- 
age (eg., galactic mask). We not only present the most general expression 
including non-uniform sky coverage as well as a non-circular beam that can 
be numerically evaluated but also provide elegant analytic results in interest- 
ing limits. As CMB experiments strive to measure the angular power spectrum 
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with increasing accuracy and resolution, the work provides a stepping stone to 
address a rather complicated systematic effect of non-circular beam functions. 
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